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ABSTRACT 

The unidentified X-ray source RX J0720.4-3125 is a candidate isolated 
neutron star showing evidence for pulsed emission with an 8.39 s period and a 
spectrum consistent with a blackbody at kT = 80 eV (Haberl et al. 1996, 1997). 
We show that this source is most likely an old isolated neutron star accreting 
from surrounding media. We then argue that unless it was born with a long 
spin period {Pi ^ 0.5 seconds) and weak field (R ^ 10^°G), the magnetic field 
on this star must have decayed. With R ~ 10^^ G, we find decay timescales 
10'^ yrs for power law decay or 10® yrs for exponential decay. A measured 
period derivative P ^ 10“^® ss“^ would be consistent with an old accreting 
isolated neutron star. Both power law and exponential decay models can give a 
P ~ 10“^® ss“^, though a P substantially less than this would be indicative of 
exponential field decay. 

Subject Headings: pulsars: general — stars: individual (RX J0720.4-3125) — 


stars: magnetic fields — stars: neutron — X-rays: stars 



1. Introduction 


There are expected to be ~ 10®-10® isolated neutron stars created in the Galaxy 
since its formation (e.g., Shapiro & Teukolsky 1983). The vast majority of these stars are 
expected to have spun down from their initial short spin periods and to have long ceased 
being active pulsars. If, however, they can accrete from surrounding media, they can 
become visible as sources of soft quasi-thermal X-rays (Shvartsman 1970; Ostriker, Rees, 
Silk 1970; Treves & Colpi 1991; Blaes & Madau 1993; Madau & Blaes 1994). Compared to 
the X-ray luminosity, these objects are expected to be very weak in the optical and infrared 
(i.e., typically Lx ~ Ltot ~ 10^° erg/s and Lopt,iR/Lx 1) regardless of whether this 
emission originates from the stellar surface or from the surrounding photoionized nebula 
(cf. Blaes & Madau 1993; Blaes et ah 1995). A detection of such old neutron stars would 
signihcantly advance our understanding of their spin and magnetic field evolution. The 
magnetic held evolution of isolated neutron stars is a major unresolved issue in compact 
object astrophysics. Theoretical studies lead, on the one hand, to exponential (Ostriker & 
Gunn 1969) or power law (Sang & Chanmugam 1987; Goldreich & Reisenegger 1992; Urpin 
et ah 1994) forms of held decay, to little or no decay within the age of the universe (Romani 
1990; Srinivasan et ah 1990; Goldreich & Reisenegger 1992) on the other. Statistical 
studies based upon the observed isolated radio pulsars give equally equivocal results (Lyne, 
Anderson, Salter 1982; Narayan & Ostriker 1990; Sang & Chanmugam 1990; Bhattacharya 
et al. 1992; Verbunt 1994) owing in large part to the difficulty in treating strong selection 
ehects (e.g.. Lamb 1992). 

In this paper, we hrst argue that the soft X-ray source RX J0720.4-3125 is an isolated 
neutron star accreting from its surroundings. Unlike the two other candidate old neutron 
stars (cf. Stocke et al. 1995; Walter et al. 1996), this source shows evidence for an 8.39 s 
rotation period — longer than in any known radio or y-ray pulsar. From this, we argue that 



4 


the magnetic field on this star must decay on timescales 10^ yrs if it was born spinning 
rapidly. A measurement of the period derivative P would help test this model. 


2. Observations 

The source RX J0720.4-3125 is an unidentified soft X-ray source seen by the Einstein 
IPC (Image Proportional Counter), EXOSAT LE (Low Energy Detector), and, most 
recently, by the ROSAT PSPC (Position Sensitive Proportional Counter) and HRI (High 
Resolution Imager). The data on this source is as follows (Haberl et al. 1996, 1997); The 
ROSAT PSPC count rate (0.1 - 2.4 keV) is 1.6 cts/s, and based upon earlier detections by 
Einstein and EXOSAT, is steady (no more than ±10% variation) over many years. The 
X-ray spectrum is best fit by a black body of kT = 80 eV with a hydrogen absorption 
column density of Nh = 1.3 x 10^*^ cm“^. With this spectrum, the count rate corresponds 
to an unabsorbed photon energy flux of i%(0.1 — 2.4 keV) ~ 1.7 x 10“^^ erg/cm^/s (K. 
Arnaud, priv. comm.). The source luminosity is then 


Lx = L{0.1 — 2.4 keV) = 1.9 x lO^^d^QQ erg/s, (1) 

where d = lOOdioo PC is the distance to the source. 

In all pointed ROSAT observations, there is a periodic modulation in the X-ray flux 
with an 8.39 s period (Haberl et al. 1996, 1997). We interpret this as the rotation period 
of the source. Since the ROSAT observations are separated by ~ 3 years, this indicates 
that the pulsed emission is steady. Two pointed ROSAT HRI observations (Haberl et 
al. 1996, 1997) secured the position of the source to be (J2000) a = 7h, 20m, 24.90s; 

S =-31° 25' 51.3" (with ±3" uncertainty). The corresponding Galactic coordinates are 
/ = 244°, b = —8°. Optical observations at the South African Astronomical Observatory 
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failed to detect an optical counterpart down to a limiting magnitude of V ~ 20.7, thereby 
placing a lower limit on the X-ray to optical flux ratio of ~ 500 (Haberl et ah 1997). 

The observational evidence points consistently to an isolated neutron star as the source. 


3. The distance to RX J0720.4-3125 

We estimate the distance to the source from the (low) hydrogen column density 
(Nff = 1.3 X 10^° cm“^). For the first 150 pc, the line of sight to this source cuts 
through the Local Bubble where the mean hydrogen density is uh ~ 0.05 cm“^ (Welsh et 
ah 1994, their Figure 3). Beyond the Local Bubble, the mean hydrogen density increases 
substantially to hh ~ 0.5 cm“^ near the Galactic plane (Dickey and Lockman 1990). Taking 
an empty Local Bubble gives a rough upper bound to the distance of about 250 pc. 

Given the very nonuniform matter distribution in the local interstellar medium (e.g., 
Welsh et ah 1994), the actual distance could be much less than 250 pc. For instance, if the 
source intercepts diffuse cirrus, its distance could be closer to ~ 100 pc (cf. Wang and Yu 
1995). A very conservative but strict upper bound on the distance is set by requiring the 
hot spot area to be much less than the star’s surface area for pulsations to be observed. 
This implies dioo 5.3i?6, where R = lO^Re cm is the stellar radius. For dehniteness, we 
adopt 100 pc throughout this work as the source distance and scale our results to this value. 
Our conclusions regarding magnetic held decay are not sensitive to the distance estimate. 


4. Arguments against a young neutron star/pulsar 


An active pulsar’s spin-down power is (e.g., Shapiro & Teukolsky 1983) Er = 
where B is the dipole magnetic held strength at the polar cap. (We took sin a = 1, 


3c3p4 , 

where a 
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is the angle between the rotation axis and magnetic dipole moment [cf. Goldreich & Julian 
1969; Verbunt 1994].) 


For young radio pulsars (age~ lO'^-lO® yrs), 3> Lx (cf. eqn [H; e.g., 


Ogelman Sz Finley 1993), yielding B 12 3> IdOdioo-Re ^-^ 8 . 39 , where P = 8 . 39 P 8.39 
s is the current observed pulsar period and B = 10^^ P 12 G. This qualihes the 
source as a “magnetar” (Duncan & Thompson 1992). The spin-down rate is 
P (-ss“^) = 2.4 X 10-16 ^ 12-^6 6 X 10“^^ dfoo-Plsg-^ 45 ^) where / = 10^^145 g- 

cm^ is the star’s moment of inertia. The consequent young spin down age — 

T'sp = ^ "C 2.2 X 10^ 6^^0-^8739-^45 FI'S — and close proximity (~ 100 pc) argues against 
a young neutron star/pulsar origin. This object lies well above the (extrapolated) 
observed radio pulsar death line (P 12 > e.g., Chanmugam 1992). Adopting Lradio 

(mJy-kpc^) = 4 X 10® P^^^/P (cf. Narayan & Ostriker 1990) yields a radio luminosity at 400 
MHz of L 400 ~ 4300 mJy. Unless our line-of-sight falls outside its radio beam 

(J. Cordes, priv. comm.), then even given the large spread (up to a factor of 100) in the 
actual radio luminosities about the best-£t Lradio{P, P) (e.g.. Lamb 1992), such a bright and 
persistent radio pulsar should have already been detected (cf. Taylor, Manchester, Lyne 
1993). The absence of radio emission further argues against a young neutron star/pulsar 
origin. 


If the source is a Geminga-type ( 7 -ray loud, radio-quiet) pulsar, then ^ Lx- 

Adopting L-y/Lx ~ 10^ as for Geminga (Halpern & Holt 1992; Swanenburg et ah 1981; 
Thompson et al. 1977) yields Tgp ~ 200 yrs, which argues against a Geminga-type pulsar. 
This is consistent with lack of detection by EGRET instruments aboard the Gompton 
Gamma-Ray Observatory (Haberl et al. 1997). 
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5. Accreting old neutron star and field decay 

If the source is an isolated neutron star accreting from the surrounding medium, its 
mass accretion rate will be governed by the rate at which material is captured within the 
star’s gravitational radius = 9.3 x 10^^ Mi. 4147 ^ cm, and is given by (cf. Bondi 

1952) Mil = facci where M = lO^^Mn g/s, M = 1.4Mi.4Mq is the stellar 

mass, riH is the hydrogen number density of the medium (assumed to have solar abundance), 
V = + = 2 OV 20 km/s with v being the star’s speed relative to the ambient medium, 

Cs being the sound speed in the ambient medium, and face is a factor that accounts for the 
microphysics of the accretion process. If the accretion is adiabatic, face ~ 1 (Bondi-Hoyle 
accretion). If preheating of the incident flow is important, face could be much less than 
unity (Shvartsman 1970; Ostriker et al. 1976; Blaes, et ah 1995; Wang & Sutherland 1997). 
Combining eqn (|^) for Ltot and Ltot = GMM/R gives V 20 = 1-1/i/c Taking 

Cs ~ 10 km/s then gives v ^ 2cs, that is, a slowly moving neutron star. In this case, 
accretion occurs in a quasi-spherical manner near the star (i.e., for S> r S> Alfven radius 
[see below]; cf. Hunt 1971). The condition n > 0 also gives another distance constraint; 
dioo < ^■^facc'^]i‘^i with the upper limit corresponding to n = 0. 

For material to reach the surface of the rotating magnetized star, the star must have 
spun down sufficiently so that accreting material can overcome the centrifugal barrier 
(cf. Illarionov & Sunyaev 1975). The neutron star’s spin evolution is divided into three 
phases (see, e.g., Blaes & Madau 1993; Lipunov, Postnov, & Prokhorov 1997). In the 
first (dipole) phase, the star is an active pulsar and spins down by magnetic dipole 
radiation; —Cl = —Cldip = = 6.2 x 10“^®i?g/^^Hi2(t)^H(f)^ (s“^), where fl = 2n/P 

is the star’s angular velocity. This phase ends when the ram pressure of the ambient 
material (~ pV"^) overcomes the pulsar wind pressure (~ Er/{ cA-nr"^)) at ~ Vg so that 
matter can now enter the star’s magnetosphere. This happens when the star’s period 




P > Pq = 4:AMi l^‘^ b12^, s. In this second (propeller) phase, material enters 
the corotating magnetosphere and is stopped at ~ r^, the Alfvenic radius, where the energy 
density in the accretion flow balances the local magnetic pressure. This radius is given by 
ta = 1.5 X 1 ^%^^ cm, where /i = BR^/2 = 10^°/i3o G-cm^ is the 

neutron star’s magnetic moment. 

Further penetration cannot occur owing to the centrifugal barrier, that is, ta > Tco (cf. 

-1 /o 

Illarionov & Sunyaev 1975), where Vco = (^) = 6.9 x lO^M^^^Pg gg cm is the corotation 

radius. Since the accretion is quasi-spherical, material falls in initially with negative total 
energy (material is bound) and roughly zero angular momentum. 

Spin down in this phase occurs via propeller and magnetic dipole spin-down; 
fl = flprop + ^dip- If the infalling material cools efficiently and attaches itself 
to field lines at around r^, the star expels the material once it spins up the 
material to the local escape velocity at ~ ta- By angular momentum conservation, 
-SiU ~ = 2.5 X (s-"). If the 

material does not cool and/or thread the field lines efficiently, the hot gas will remain bound 
in a “cocoon” at ~ ta- As the underlying magnetosphere rotates supersonically shearing 
through this “cocoon,” the consequent shock heating expels material from the star (Illarionov 
& Sunyaev 1975). By energy conservation, -Vt^rop ~ 7 ^ = ~ 

where the final expression is typical for a star just entering the propeller phase. In this 
limit, only a fraction of the energy given to escaping material goes into azimuthal motion, 
so —(iprop ^ —^prop- We therefore take tlprop = fp^'prop: where 0.01 ^ /p ^ 1 is taken as 
a freely adjustable parameter. The exact value of fp requires time-dependent numerical 
simulations (e.g., Wang & Robertson 1985). If fp ~ 1, Clp^op 3> Cl^ip throughout this phase. 

Propeller action continues until ta < Tco, when the centrifugal barrier is removed and 
polar cap accretion ensues (e.g.. Lamb et ah 1973; Davidson & Ostriker 1973; Arons & Lea 
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1980). This occurs when P > Pa = s- In quasi-spherical 

accretion, no net torque is exerted on the star (due to accretion). The mass loading of the held 
lines, however, now becomes important in spinning down the star by increasing the moment of 
inertia of the star + corotating magnetosphere systemic (see, e.g., Mestel 1990 for a discussion 
of this effect in normal stars). For ta ^ R-, angular momentum conservation (^(Jf2) = 0) 
gives -tlkrk = ~ = 4.5 X (s'^). (In 

the limit of a weakly magnetized star where ta < R, ta ^ R in the above expression 
and —(ibrk —*■ % fl-) The total spin-down rate in this hnal (accretion) phase is then 
n = flbrk + ^dip, although, quite generally, flbrk ^ ^dip- 

For accretion to occur, we require ta < Tco, that is, P > Pa. For pulsed emission, that 
is, polar cap accretion, we require ta 3> R. Combining these, the present day surface held 
must satisfy 


« Bn < B„„,n = 9.3 x lO-U/'f 


- 3/2 1/2 j-,7/6 
^ 8 . 39 - 


( 2 ) 


This condition rules out the possibility of measuring the held directly through cyclotron 
(emission) line observations, which require B 12 ~ 1 (Nelson et ah 1995). 

Assume hrst that the star’s held does not decay, that is, its initial held Bi is the same 
as the current held and satishes eqn (^. We obtain an upper bound to the period Pq by 
using B < Bcrit (cf. eqn [|]); 


p„ < =o.42Mf.y/yv2i'/‘‘ s. 


( 3 ) 


thank Eve Ostriker for pointing out the potential importance of this ehect. 
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If the star was born with Pi < Pq < Peru-, h must hrst dipole spin down to Pq. This 
takes tdipfl = 6.3 x lO'^ P^B^^Re^hb yi’s > 1.3 x 10^^ h 5 R 8 l 9 ^''^H^^ 2 d‘^ where we 

have used B < Beru (cf. eqn p|) to arrive at the inequality. The star thus spends longer 
than a Hubble time (~ 10^° yrs) just in the hrst (dipole) phase. 


If the star was born with Pq < Pi < Peru-, h goes directly 
to the second (propeller) phase. Propeller spin down to Pa takes 
tprop,a 1-1 X 10 Jp JPJ-IA Rq acc ^20 -^12 J ^ 

10^° f-^ yrs, where we have used B < B^ru (cf. eqn 

0]), Pi < Perit (cf. eqn [^), and Pi <C Pa to arrive at the inequality. The star thus spends 
longer than a Hubble time in this phase. (If /p -C 1, magnetic dipole spin down will be 
important initially, but propeller spin down dominates eventually.) 


We conclude that if the star was born with Pi ^ Pent = fUc^ 2 o^^'^ s, then to 

enable sufficient spin down to allow accretion onto the star, it must have been born with a 
stronger held than at present, that is, the stellar magnetic field must have decayed. 


Of course, if the neutron star is born slowly rotating, i.e, with Pi 3> Peru, then the 
star can spin down within a Hubble time to 8.39 s without requiring magnetic held decay. 
However, it is unclear how neutron stars can form with such long initial spin periods and 
weak magnetic helds (cf. eqn P|). One possibility is that RX J0720.4-3125 actually evolved 
from a high mass X-ray binary system (cf. Haberl et ah 1997). 0 We argue here, however, 
that conventional isolated neutron star models can account quite well for the properties of 
RX J0720.4-3125. 


In general, the star’s spin and magnetic held history is determined by Pj, Pj, the held 


^It is also possible that no accretion is involved and that RX J0720.4-3125 is powered 
instead by some internal heat source (see, e.g., Thompson & Duncan 1996). 
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decay law, and decay timescale td- The dependence on Pj, however, is very weak whenever 
Pi -C Pa, which is generally believed to be the case (e.g., Narayan & Ostriker 1990). 

In Fignre 1, we illnstrate sample evolntionary tracks in B-P space for two magnetic 
held decay laws. For the solid cnrve, the stellar held is assnmed to decay as a power law, 
i.e., B{t) = Bi/{1 + t/td) (e.g., Narayan & Ostriker 1990; Sang & Chanmngam 1987), 
with td = 3.8 X 10^ yrs. For the dot-dashed cnrve, exponential decay is assnmed, i.e., 
B{t) = Bi exp{—t/td) (e.g., Ostriker & Gnnn 1969), with = 4 x 10® yrs. To constrnct 
these tracks, an isolated nentron star with M = 1.4Mq and P = 10 km is assnmed to be 
born with Pj = 0.01 s < Peru, Bi = 10^^ G, and is assnmed to be moving at 20 km/s 
throngh a medinm with uh = 1 cm“® and Cs = 10 km/s. Spin down spans the magnetic 
dipole, propeller, and accretion phases. For the propeller spin down rate, we took fp = l 
and for the accretion rate, we took face = 1- In the power law decay model, the star enters 
the accretion phase at 5.7 x 10® yrs and spins down to 8.39 s in 6.2 x 10® yrs. For the 
exponential decay model, these nnmbers are 1.9 x 10® yrs and 2.1 x 10® yrs, respectively. 

For Pj_i 2 ~ 1, we hnd td ^ 10^ yrs for power law decay models (cf. Urpin et ah 
1994), while td ^ 10® yrs for exponential decay models. These results are consistent with 
several recent analyses of pulsar statistics and held decay (e.g., Wakatsuki et al. 1992; 
Bhattacharya et al. 1992; Lamb 1992; Urpin et al. 1994; Verbunt 1994), and with the 
analysis assuming power law decay in Narayan & Ostriker (1990). For Bi 12 ^ 0.1, the star 
cannot enter the accretion phase within a Hubble time even for the most favorable case of 
a stationary (n = 0) star. If /p -C 1, the star will not be able to enter the accretion phase 
and spin down to 8.39 s in a Hubble time unless n <C 20 km/s (cf. Blaes Sz Madau 1993). 

From flbrk and B < Beru (cf. eqn [Q), we obtain an upper bound to the current spin 


down rate; 


12 


P 


< 1.8 X 


r5/3 D— 2 


NOW 


-3 




/< 


p7/3 

acc-rg,39 


SS 


-l^ 


(4) 


In both models shown in Figure 1, the 8.39 s period is reached shortly after the star enters 
the accretion phase, and P\jqow ~ 10“^® ss“^. 

It is evident from Figure 1 that P in the exponential decay model asymptotes after 
entering the accretion phase while P for the power law decay model continues to increase 
monotonically. This is because in exponential decay models, the field, and hence ta (lever 
arm) decreases more rapidly at late times, whereas the field in power law decay models 
persists longer. If neutron star helds undergo exponential decay and the 8.39 s period in 
RX J0720.4-3125 is the star’s asymptotic period, then we expect P\now ^ 10“^® ss“h 

We have so far assumed that the stellar magnetic held decays indehnitely. We 
emphasize, however, that once the stellar held has decayed sufficiently to enable accretion 
(cf. eqn [2]), no further decay is required by the observations. Thus, for instance, the 
decaying held may level out at late times to a steady hnite long-lived value (e.g., Kulkarni 
1986; Romani 1990). 


6. Conclusions 

We argue in this paper that RX J0720.4-3125 is an old isolated neutron star, situated 
at ~ 100 pc, that has spun down long past the active pulsar stage and is now accreting 
matter from the interstellar medium. Unless the star was born with an unusually long 
period (Pj 0.5 s; cf. eqn |^) and weak magnetic held (Pj ^ 10^° G; cf. eqn 0), such spin 
down is only possible if the star’s held at birth was much stronger than at present, that is, 
the star’s held must have decayed. Our analysis gives long decay timescales; 10^ yrs for 
power law decay and 10® yrs for exponential decay, assuming Pj ~ 10^^ G. 
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A P ^ 10 ss ^ would be consistent with an old accreting isolated neutron star. Both 
power law and exponential decay models can give a P 10-16 ss-i. A P < 10-16 ss-i, 
however, would be indicative of exponential held decay. In addition, an accreting isolated 
neutron star should be surrounded by an extended photoionized nebula (cf. Shvartsman 
1971; Blaes & Madau 1993; Blaes et al. 1995). This low surface brightness nebula should 
be dominated by Ha in a compact inner zone (~ lOi^ cm; e.g., Blaes et al. 1995) and be 
rich in mid to far infrared metal forbidden lines (e.g., [NeII]12.8/i, [SiII]35/i) in an extended 
outer zone (~ lOi® cm). The results of deep spectroscopy can therefore be used as an 
additional test of the model presented here. 

We thank Gunther Basinger for pointing out this source to us, and we thank Doug 
Hamilton, Eve Ostriker, Jim Stone, Mark Wolhre, and Sylvain Veilleux for valuable 
discussions and a critical reading of the manuscript. We also thank an anonymous referee 
for insightful comments and suggestions. This work was supported in part by NASA 
Astrophysics Theory Program grant NAG5-3836. 
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Fig. 1.— The spin period and magnetic field evolutionary track for a neutron star 
(M = 1.4 Mq, i? = 10 km) born with period Pi = 0.01 s and surface dipole field strength 
Bi = 10^^ G. The star is assumed to be moving at n = 20 km/s through a medium (solar 
abundance) with density uh = 1 cm“^ and c* = 10 km/s. The star’s field undergoes (solid 
curve) power law decay {B{t) = 5j/(l Tt/t^)) with td = 3.8 x 10^ yrs, or (dot-dashed curve) 
exponential decay {B{t) = Biexp{—t/td)) with = 4 x 10® yrs. Dashed line gives the 
observed radio pulsar death line {Bu = 0.2P^). For power law decay, the star drops below 
the death line at t ~ 3.5 x 10® yrs after its birth when P = 0.73 s. Propeller spindown 
begins at t ~ 1.3 x 10® yrs when P = Pq = 0.76 s. The star enters the accretion phase at 
t ~ 5.7 X 10® yrs when P = 6.96 s. Spin down to the observed 8.39 s period occurs after 
6.2 X 10® yrs. For exponential decay, the corresponding times and periods are (3.1 x 10® yrs, 
1.56 s), (7.5 X 10® yrs, 1.74 s), (1.9 x 10® yrs, 7.86 s), and (2.1 x 10® yrs, 8.39 s), respectively. 
Each of these events is marked with an open circle on the curve for the power law decay 
model (solid) and an open triangle on the curve for the exponential decay model (dot-dash). 





